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A recently discovered material, 6H-B-Ba3NiSt>209 was found to display unusual low-temperature phe- 
nomenology, interpreted as a quantum spin liquid with spin S = 1 on a triangular lattice. We study a spin 5 = 1 
exchange model on an AB stacked triangular lattice near its quantum paramagnet-to-spiral transition, driven by 
easy-plane single-ion anisotropy. We demonstrate that the frustrated inter- and intra-layer exchanges induce 
contour lines of low-energy excitations that lead to a broad crossover regime of linear-temperature dependence 
of the specific heat. Based on this and various other predictions, we argue that the observed phenomenology can 
be understood in terms of a conventional picture of a proximity to this frustrated critical point. 

PACS numbers: 71.70.Ej,71.70.Gm,75.I().-b 



Quantum spin liquids (QSLs) are Mott insulators that re- 
main magnetically disordered down to zero temperature, and, 
as we use term here, are exotic states of matter characterized 
by properties such as quantum number fractionalization, topo- 
logical order, and gapless excitations in the absence of spon- 
taneously broken symmetry. The realization of QSLs in theo- 
retical models has been well established! 1|, an d a number of 
materials have emerged as promising candidates|2]-(9]. How- 
ever, there is no direct confirmation of QSL in any of these 
systems, and alternative explanations now exist for some QSL 
candidates Ifl0ifl4l . 

Many QSL candidates share a rough phenomenology: they 
are electrical insulators, but with thermodynamic properties 
similar to those of a metal. In particular, many of these sys- 
tems have a constant low-temperature spin susceptibility, and 
a linear-temperature dependence of the low-temperature spe- 
cific heat. Theoretical attempts to explain this behavior usu- 
ally invoke spin-^ fermionic spinons with a constant density 
of states (DOS). In this Letter, we propose the first (to our 
knowledge) alternative explanation for this phenomenology 
that does not invoke substantial quenched disorder. 

Recently the compound 6H-B-Ba 3 NiSb 2 9 (6H-B) has 
been proposed as a QSL candidate !!?! . This system has mag- 
netic ions Ni 2+ forming triangular layers with spin-1 local 
moments. The Curie-Weiss temperature is — 75.5K and no 
sign of magnetic ordering is detected down to 0.35K, indicat- 
ing a strong frustration. The system exhibits the QSL phe- 
nomenology described above, with a linear-T specific heat 
and constant spin susceptibility at low temperature s |[T5l . To 
account for the experiments, Refs.fT^landfTTlproposed QSLs 
with fermionic spinons. In contrast to these interesting pro- 
posals, in this Letter we argue that the 6H-B data can be un- 
derstood without invoking QSL physics. We propose that the 
putative QSL behavior arises as a crossover tied to the prox- 
imity of a quantum critical point (QCP) between spin spirals 
favored by the frustrated exchange, and a quantum paramag- 
netic (QP) phase, favored by a single-ion anisotropy (SIA). 

More specifically, in a mean-field treatment we find the 
dispersion of spin excitations has the schematic form = 
V '/i(fc)/2(k)- While generically, including at the QCP, there 
is no special relationship between the functions fx and / 2 , in 



a broad parameter regime near the QCP fx and fi are approx- 
imately proportional. This leads to an enhanced DOS, and, 
due to the presence of a degenerate contour of low-energy ex- 
citations, a linear intermediate-temperature specific heat. This 
behavior follows from the form of the dispersion, and is ex- 
pected to be robust beyond mean-field theory (MFT). The 
microscopic ingredients for this behavior are SIA combined 
with comparable Ising and transverse antiferromagnetic ex- 
change. Therefore, we expect that such a deviation from a 
generic dispersion, accompanied by anomalous intermediate- 
temperature thermodynamics, should be common in S > 1/2 
antiferromagnets where the crystal structure admits a SIA. 
More broadly, there are certainly many mechanisms by which 
generic behavior may be pushed down to very low tempera- 
tures, and the resulting regimes of anomalous intermediate- 
temperature behavior may be important in various situations, 
perhaps even in other QSL candidates. 

In 6H-B, the Ni triangular layers have an A-B stacking with 
the lattice sites on one layer projecting to the centers of the tri- 
angle plaquettes on the two neighboring layers (Fig. [TJ. Our 
minimal model includes the interlayer and intralayer spin ex- 
change and a SIA. Treating the two neighboring triangular 
layers as the two sublattices of a honeycomb lattice, we view 
the system as a multilayer honeycomb lattice (Fig. [TJ. There- 
fore, when the exchange is dominant and frustrated, the clas- 
sical ground state is highly degenerate. Quantum fluctuations 
lift the degeneracy and favor coplanar spiral orders. A strong 
easy-plane SIA favors a QP state, which is separated from the 
ordered state by a QCP. We propose that 6H-B is close to this 
QCP, and may lie either on the QP or magnetically ordered 
side. The constant spin susceptibility arises from the explicit 
breaking of spin rotational symmetry by the SIA, and the pow- 
der nature of the samples. More notably, we interpret the ob- 
served broad linear-T specific heat in terms of the dispersion 
for spin excitations, as discussed above. 

Model — Although the interlayer exchange path goes 
through one more oxygen than the intralayer coupling, the 
multiplicity of the former path is larger than the latter. More- 
over, in a structurally similar material 6H-A-Ba3NiSb20g 
with long-range magnetic order, the magnetic specific heat 
at low temperatures is observed to behave as C V (T) ~ T 3 , 
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FIG. 1. (Color online) Two adjacent triangular layers of 6H-B (left) 
can be viewed as a single-layer honeycomb lattice (right). Ji (J2) 
is the interlayer (intralayer) exchange| 19]. Ji is also the exchange 
between adjacent honeycomb layers. 



which indicates a non-negligible interlayer coupling. That 
may thus be important in understanding the properties of 6H- 
B, but is not required for our theory. The resulting exchange 
is given on the triangular multilayers by the Hamiltonian, 
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The first sum is for interlayer exchange between nearest- 
neighbor (NN) sites on neighboring A and B layers, and 
the second sum is for intralayer exchange between NN sites 
within the same layer. As illustrated in Fig. [T] the interlayer 
(intralayer) exchange on a triangular bilayer can be viewed 
as the nearest-neighbor (next-nearest-neighbor) exchange on 
a single honeycomb layer. In contrast to Ref. [161 we do 
not include the biquadratic exchange, which we expect to be 
strongly subdominant to H ex . 

The space-group symmetry P63II1C of 6H-B restricts the 
SIA to be "Hani = D X^(<S?) 2 witn z-axis normal to the tri- 
angular layers. Since an easy-axis anisotropy is more likely 
to favor magnetic order, so we expect easy-plane anisotropy 
(D > 0) for 6H-B, where such order is not observed. 

Our model thus contains two competing terms, exchange 
and SIA, T~L = W ex + H dn \. Implementing high-T series 
expansion, we extract the Curie- Weiss temperature, finding 
that 9g w = -D/3 - 4 J and 9^ w = D/Q - 4 J (where 
J = Ji + J 2 ) for field applied along and perpendicular to the 
z axis, respectively. With a powder sample in experiment l[T5l . 
a powder average gives Qq W = — 4 J that is independent of D. 
Furthermore, with Weiss-MFT we demonstrate that saturation 
temperature of spin susceptibility (observed to be ~ 25K| 15 1) 
is set by D, that is thus comparable to J. 

For the Hamiltonian H, when the SIA dominates with 
D J, the ground state is a uniform QP state with \S Z =0) 
at each site. In the opposite limit of dominant exchange, we 
expect the ground state to be magnetically ordered. Luttinger- 
Tisza method[20| gives the classical ground state spin con- 
figurations with the ordering wavevector q z = and spins 
lying in the xy plane. When Ji > 3 J2, the classical ground 
state is a usual Neel state. When J\ < 3J 2 , the classical 
ground state is degenerate with degenerate spiral wavevectors 
q ± = (q x ,q y ) satisfying J2{b} cos(q_i_ • b) = (^) 2 - 3, 
where {b} are 6 next-nearest-neighbor lattice vectors of the 
honeycomb lattice. The degenerate wavevectors form contour 



curves in momentum space. Moreover, with vanishing J\, this 
spiral reduces to the familiar 120° state of decoupled triangu- 
lar layers. Quantum fluctuations lift the degeneracy of these 
classical spin spirals, selecting states characterized by a dis- 
crete set of q's around which the quantum zero-point energy 
is minimized. The spiral ground states favored by the quan- 
tum fluctuations do not vary upon introducing the SIA. The 
optimal spiral wavevectors are given by 1 18 1 
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and their symmetry equivalents. 

Starting from the magnetically ordered phase, the existence 
and properties of the phase transition can be analyzed within a 
Weiss-MFT. We decouple the exchange into an effective Zee- 
man field which is then self-consistently determined for each 
sublattice. We parameterize the spin order as, 

S A (r) = A/[cos(q • r)x + sin(q ■ r)y], (4) 
S B (r) = A/[cos(q • r + 6)x + sin(q • r + 0)y], (5) 

where 8 is the relative phase between two sublattices and 
depends on j 1 , and M is the order parameter to be deter- 
mined. This parameterization describes both the Neel state 
for J 1 > 3 J 2 and the spin spirals for J\ < 3 J 2, with the 120° 
state as the limiting case of the decoupled triangular layers. At 
zero temperature MFT yields that in the vicinity of the QCP 
the order parameter is M = v2(l — with the critical 

anisotropy parameter D c = 12(Ji — J2) for the Neel state 
when J 1 > 3 J 2, D c = 6 J 2 for the 120° state at vanishing J l5 

2J 2 

and D c = 6 J 2 + -j^- for the spin spirals when J\ < 3J 2 - We 
expect that as usual Weiss-MFT overestimates D c (Fig. [2]). 

Within Weiss-MFT, in the QP phase, the zero-temperature 
spin susceptibility \ z = (fields along z-axis). For fields in 
xy plane, the spin susceptibility saturates to a constant \o — 
2 ^D+i2j • T ne P ow der average gives the zero-temperature 
susceptibility — %Xo /3- 

Rotor MFT — It is convenient to model this easy-plane sys- 
tem with rotor variables, by introducing an integer-valued 
field m and 27r-periodic phase variable (pi, which satisfy 



bi, rij\ — iSij. With the mapping (Sf 
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the rotor Hamiltonian reads 
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where takes J\ (J2) for NN interlayer (intralayer) bonds. 
Although rii only takes the values of ±1, in the spin model, 
due to the substantial anisotropy D, we expect that relaxing 
this restriction is unlikely to have significant effects. 

Using the coherent-state path integral, we integrate out the 
field rii and obtain the partition function, 



Z = /2?$Z>Ae-' s - < S,/«h-A*(l*il"-i) 
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FIG. 2. Zero-temperature phase diagram determined from the SPE. 
Dashed curve indicates the location where D = Ao/2, which is 
important in the discussion of T-linear C v (T) below. 



where S = / dr£ k (4DI + 2J k )-X^A$ ! , ) -k + 
J2ij Jij^i^j Wltn = e**'. .7k is the 2 x 2 exchange cou- 
pling matrix in momentum space, /x, are the sublattice in- 
dices, and I is a 2 x 2 identity matrix. The constraint |$;| = 1 
is enforced by the Lagrange multiplier A;. We proceed by 
a saddle-point approximation. Assuming i\ = (3A(T) at 
the saddle point, we integrate out the $ field and obtain the 
saddle-point equation (SPE) for A(T) in paramagnetic phase, 
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where ubz = is the Brillouin zone (BZ) volume, s+ k 

v3 ' 

J2 E {b} cos(k-b)±2| Ji cos(^)|^3 + S {b} cos(k ■ b) 
the eigenvalues of J7k, and e±^ are the two spin excitations 
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When the left-hand side of the SPE is less than 2 for any 
A(T), the rotor is condensed which signals magnetic order- 
ing. Therefore, besides the transition temperature from the 
high-temperature paramagnetic phase to the low-temperature 
spin spirals, we also obtain the critical D c that separates spin 
spirals from QP phase and the zero-temperature phase dia- 
gram(Fig. [2]). As expected, D c obtained here is smaller than 
the one determined previously from the Weiss-MFT. In partic- 
ular, D c / J is minimal at Ji = J 2 corresponding to the largest 
frustration at this point. Right at the QCP and zero tempera- 
ture, A(0) = Ao = 3J2 + J1/J2 an d the low-energy mode 
e_.k develops gapless excitations. As shown in Fig. [3j the 
momenta of the gapless excitations form contour lines that are 
identical to the ones of degenerate classical ground state spiral 
wavevectors. Moreover, as J\/J% increases from 0, the con- 
tour lines around the BZ comers gradually expand and meet 
at M when J\ = J 2 ■ 

Near the QCP with T < J — A(T) increases with T and 
we define A(T) = A + A X (T). The excitation e±(k) 
picks up a self-energy via the T-dependence of A(T). By 
numerically solving the SPE, we find that, near the QCP 




FIG. 3. (Color online) The evolution of the low-energy excitations 
in k x -ky plane with k z — at the QCP. The parameters used in the 
figures are (a) Ji = 1.5J 2 ,D C = 1.36 J, (6) Ji = J 2 , D c = 1.23 J, 
(c) Ji = 0.8J 2 ,D C = 1.28J, (d) Ji = 0,D C = 2.01J. The low- 
energy gapless contours are marked with bold black lines in (a-c), 
while in (d) the low-energy gapless points are marked with black 
dots. Lattice constants are set to 1. (e) is the BZ of a honeycomb 
lattice. For Ji > J2, the contour line is centered in the middle of BZ. 
For J\ < J2, the contour lines are centered around and eventually 
shrink to the corners of BZ in the limit Ji — > 0. The "O" in (6) 
correspond to M in (e). 



Ai(T) oc T 2 for T < J. This is also supported by an an- 
alytical argument (Supplemental Material |fT9ll ), and holds in 
the quasi-2D limit J\ <C J%. This immediately leads to the 
internal energy E oc T 3 and hence C v oc T 2 , for T <C 
with n an energy cutoff. This low -temperature T 2 -C v regime 
is confirmed numerically in Fig.Qa). 
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FIG. 4. (Color online) The C v -T plots in the paramagnetic phase. 
In (a), from top to bottom, J\ — J2, 0.7 J2, 1.5J2, 1.8J2, O.3J2, 
with D w D c and D = 1.24J, 1.31 J, 1.36J, 1.47J, 1.58J, 2.01J, 
respectively. In (b), J 1 = 0.5J 2 , D = 1.08J, 1.17J, 1.23 J, 
1.32J, 1.41J, 1.48J, 1.55J from top to bottom. For D = 1.08J, 
1.17J, 1.23 J, 1.32J, T c = 0.10J,0.07J,0.064J,0.04J, respec- 
tively. For Ji = O.5J2, Dc — 1.41J. The dashed lines are the 
linear fits for a range of data points. Energy is in units of J. 

In Fig.Qa), we also find that, as J\ / J2 moves to the point 
J 1 = J2 from either side, the temperature range of the T 2 -C v 
regime diminishes. We attribute this to the observation that the 
zero-temperature DOS at the QCP increases with energy, then 
saturates to a roughly constant value. This saturation energy 
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(in units of J) is found to be lowest when J\ = Ji- 

Linear-T C v at intermediate T — In Fig. |4] we find an 
intermediate-temperature regime with C v rj c\T + cq. To ex- 
plain this, we first note that when D « y 2 ' , the low -energy 
excitation is approximately the square root of a perfect square: 

e_ )k w V2|s-, k + D + A(T)/2| 

2 2mj_ ifco 2m z ,fe 

If = 4p-, such dispersion on contour lines in 3D gives 
a constant DOS at low energies and, because A(T) is only 
weakly T-dependent in this case, this leads to C v oc T at low 
T. This conclusion also holds in the 2D limit J\ <C J2. As 
shown in Fig. [2] D c is slightly greater than Aa and the system 
at D = A 1 is ordered at very low T. Once the system enters 
into the paramagnetic phase, a linear-T C v is obtained (see 
the top curve in Fig.Qfe)). Moreover, this linear-T C v regime 
persists even when D is increased to or slightly beyond D c . 

Discussion — The spin susceptibility is observed to saturate 
to a constant below 25K lfT31 . which is consistent with our pre- 
diction. Experiments also find C V (T) ~ with 77 s» 1.0(1) 
for 0.35K < T < 7K. Determination of both behaviors relies 
on subtracting a magnetic impurity contribution, which has a 
significant effect below about 25K and IK for \ and C v , re- 
spectively. The subtraction procedure for C v relies on fitting 
C V (T, B = 0T) — C V {T, B = 9T) to a Schottky form appro- 
priate for a magnetic impurity contribution. This procedure 
may be unreliable, because 9T is a large field scale for 6H-B, 
corresponding in temperature units to ~ 20K, and the fitted 
difference of C v is thus expected to include a significant con- 
tribution from bulk Ni moments. While this complicates in- 
terpretation of the C v (T) for T < 1 K, we note that C v (T) is 
certainly linear for IK < T < 7K (or 0.05 J < T < 0.34 J). 

As illustrated in Fig. |4] for some parameter values, the 
range of the linear-T specific heat of the experiments is 
compatible with our results for a range of J1/J2IISI. Tak- 
ing Ji = 0.5 J 2 and D = 1.32 J (Fig. Qb)), we find the 
powder-averaged zero-temperature spin susceptibility Xo ~ 
0.0124emu/mol, which is very close to the experimental 
value 0.013emu/mol. Moreover, the cofficient 7 in C V (T) is 
found to be 204mJ/mol-K 2 while the experimental value is 
168mJ/mol-K 2 . We obtain a Wilson ratio of 4.4, not far from 
the experimental value 5.6. The agreement can be further im- 
proved by adjusting J\j J 2 and D. 

6H-B may lie either on the QP or magnetically ordered side 
of the QCP. At very low temperatures we thus expect either 



a small energy gap, or the onset of spin order. It should be 
noted that the presence of magnetic impurities may interfere 
with observation of such very-low-T behavior. To detect the 
energy gap or spiral order, NMR or /iSR measurements may 
be helpful. If 6H-B is in the QP phase, then in a single crystal 
sample we predict \ z = 0, x = const at zero temperature. 
Re f. Unconsidered a state with gapped S z — ±1 fermions and 
gapless S z = fermions forming a Fermi surface. This state 
has a spin gap, with thermodynamic properties dominated by 
the Fermi surface. The two QSLs in Ref. 17 have gapless spin 
excitations. These states are thus distinct from the QP phase of 
our proposal, which has a fully gapped spectrum. Especially 
if single crystal samples are available, inelastic neutron scat- 
tering should be able to further distinguish these proposals by 
measuring the dispersion of low-energy spin excitations; our 
prediction is depicted in Fig. [3] 

To summarize, we propose a minimal J1-J2-D model for 
6H-B and argue that its putative QSL phenomenology is 
due to the proximity to a QCP. Our theoretical prediction is 
broadly compatible with current experiments! 15]. Various fu- 
ture experimental directions are suggested. 
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SUPPLEMENTARY MATERIAL 

Here we provide a detailed analytical argument for the T 2 specific heat at low temperature near the QCP. At finite temperature, 
A(T) increases with temperature and we define A(T) = Aq + Ai (T). The spin excitation e± (k) also picks up a self-energy via 



FIG. 5. (Color online) The couplings in the multilayer triangular structure of 6H-B-Ba3NiSb20(). Ji, J2 are the interlayer and intralayer 
couplings, respectively, (a) is the three dimensional view of the lattice structure. The dashed circles on the B plane are the projected positions 
of the lattice sites from the A plane, (b) is the top view of the lattice (aslo shown in the main text). 



the temperature dependence of A(T). By numerically solving the saddle-point equation, we find that, near the QCP Ai (T) oc T 2 
at T <C J. At T <C J, the low-energy spin excitation near the contour lines can be approximated as, 



AMT)+vl ko ki+vl ko kl (11) 

where A = 4D C — 2Aq, ko is a momentum coordinate running along the contour lines, k± is normal to the tangent of the 
contour line at ko, and we have neglected the weak temperature dependence of the speeds Vj_, v 2 . Eq. ( fTTj ) is expected to be a 
good approximation for e_ (k) less than a cutoff energy CI with T <C Cl <C J. The saddle-point equation can be approximated as 

A coth( I JAA l + v 2 k 2 +ufM) 

\ 2 v - 1 - 1 z z ' + b= 2, (12) 

k ,k x ,k z 2 y /AA 1 + v 2 ± kl + v 2 k 2 

where the integral is over the region around the contour lines with |fcj_|> \k x \ ^ A, and b is the approximately T-independent 
contribution from outside this region. At low temperature, the temperature-dependent part of the integral becomes independent 
of the cutoff f2, and only depends on T via the dimensionless parameter — ^ ■ In order for the integral to be constant in 
temperature, we thus expect Ai(T) oc T 2 in the limit T <gi CI. This result immediately leads to the internal energy, which can 
be approximated as 



JAAUn+v^k 2 +v 2 k 2 , 
£~ / v , v ' - 1 - 1 z z - oc T 3 (13) 



for T ^ CI. This gives C v oc T 2 in this temperature regime. 



